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ABSTRACT 

In a previous paper we identified cores within infrared dark clouds (IRDCs). We 
regarded those without embedded sources as the least evolved, and labelled them 
starless. Here we identify the most isolated starless cores and model them using a 
three-dimensional, multi-wavelength, Monte Carlo, radiative transfer code. We derive 
the cores' physical parameters and discuss the relation between the mass, temperature, 
density, size and the surrounding interstellar radiation field (ISRF) for the cores. The 
masses of the cores were found not to correlate with their radial size or central density. 
The temperature at the surface of a core was seen to depend almost entirely on the level 
of the ISRF surrounding the core. No correlation was found between the temperature 
at the centre of a core and its local ISRF. This was seen to depend, instead, on the 
density and mass of the core. 

Key words: stars: formation - Infrared Dark Clouds 



1 INTRODUCTION 



IRDCs were initially d iscovered b y the MSX |Carev et ahl 
ll99Sl ; lEgan et al.lll998l ) and ISO jPerault et al.lll996l ) sur- 
veys as dark regions against the mid-infrared (MIR) back- 
groun d. The densest IRDCs may eventually form massive 
stars ( Kauffm ann fc Pillail bOld ). and are therefore pre- 
sumed to represent the earliest observable stage of high mass 
star formation. 

Some IRDCs contain cold, compact condensations 
called infrared dark c ores; also referred to as 'clumps' (e.g. 
Battersby et al.|[201ol ) or 'fragments' (e.g. iPeretto fc Fuller! 



20091 ). The most massive of these objects are believed 



to be th e high mass equivale nt of low-mass prestel- 

lar cores (IWard-Thompson et all Il994l ; ICarev et all |2000| ; 
iRedman etaT]|2003t iGarav et al.ll2004 ). 

Observations of IRDCs have shown them to have low 



temp eratures (T<25K; e.g. lEgan et al.ll 19981; iTevssier et ail 



20021). high den s ities (W >10 5 cm" 3 ; lEgan et all Il998l : 



Carey et alJll998t l2000h and masses betw een ~10 2 -10 4 M Q 
l|Rathborne et alj|2005fc iRagan et al.ll2009D . 

Within the IRDCs, infrared dark cores have been ob- 
served with radii ranging from 0.04 to 1.6 pc, m asses from 10 
to IOOOMq and densitie s from 10 3 to 10 7 cm" 3 (|Caray et all 
| 2004| ; lOrmel et al.ll2005l ; iRathborne et all 120051 . 120061 ; ISwiftl 
120091 ; IZhang et al.ll2011 ) . This is consistent with them being 
the precursors to hot molecular c ores, the next observab le 
stage in high mass star formation l|Rathborne et afl feOQO), 

Using CO and CS velocities, kinematic distances have 
been calculated fo r several hundred I RDCs and, by exten- 
sion, their cores (ISimon et al.l 120061 ; I Jackson et al.l [2008). 
The distribution of IRDCs in the first Galactic quadrant 
differs from those in the fourth quadrant. In the first quad- 
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Table 1. The physical properties of the modelled cores. Column 1 gives the name of the parent cloud as it appears in PF09. Columns 2 
and 3 show the Right Ascension and Declination of the cores in degrees. Column 4 is the calculated kinematic distance from the sun to 
each core (see Section l3.2l l. These distances have an uncertainty of ~20 per cent based on how accurately the CO velocity can be found . 
Column 5 shows the interstellar radiation field needed to match the model to core's SED, given in terms of multiples of the Black (1994) 
radiation field. Column 6 is the central density of the core. Column 7 is the core's semi-major axis, as measured at 250 (im. Column 8 
is the flattening radius of the core and is set to be one-tenth of the semi-major axis. Column 9 gives the asymmetry factor, a measure 
of the eccentricity of the core. Column 10 is the mass of each individual core. Columns 11 and 12 show the temperatures at the centre 
and the surface of each core, respectively, assuming a viewing angle of 8=0 ° . The temperatures and masses of the cores are calculated 
by the model (see text for more details). The x 2 values show how the goodness of fit of the model. See text for details. 



Name of the RA Dec Distance ISRF »io(H2) Umax Ro Asymmetry Mass Temperature (K) \ 2 

Parent IRDC 2000 (°) 2000 (°) (pc) (cm" 3 ) (pc) (pc) factor (M ) Centre Surface 
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3.0xl0 4 


0.9 
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208 
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-61.87 


3600 


0.7 


S.OxlO 3 


0.9 


0.09 
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3.0 


53.8 


13.5 


23.9 


1.94 
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51.4 
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322.334+0.561 
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rant IRDCs have a m ean galacto-centric distance of ~5 kpc 
(| Jackson et al.| [2008). whereas those in the fourth qu adrant 
are typically found at ^6 kpc ijPeretto fc Fullerll2010t l. This 
suggests that many IRDCs may be located in the Scutum- 
Centaurus spiral arm. The association of IRDCs to a spiral 
arm lends weight to the theory that IRDCs are the birth- 
place of high mass stars (jjackson et al.ll2008l ). 

Cores within IRDCs can be separated into different 
classes based on observational parameters. Starless cores, 
also called 'quiescent cores', are the youngest subset of in- 
frared dark cores. They are not yet undergoing any form 
of active star formatio n and so contain no MIR activity 
l|Chambers et al.ll2009h . As the cores evolve they leave the 
starless phase and begin to exhibit tracers of star forma- 
tion. First_24^xrn ejnission, and then PAH 8-/xm emission 
(see IWilcock et ail |2012| for more details). When the cores 
have reached the latter stage they ar e thought to contain a 
hyper- or ultra-comp act HII region l|Chambers et alj|2009l ; 
iBattersbv et al.ll2010t) . 

In a previous paper (IWilcock et all I2012T I. we used 
data f rom the Herschel Infra red Galactic Plane Survey (Hi- 
GAL; iMolinari et all l2010al ibh to observe the IRDCs from 
iPeretto fc Fuller! (|2009h - hereafter PF09 - at FIR wave- 
lengths. We identified 1205 IRDCs that were simultaneously 
Spitzer-dark in the MIR and Herschel-bright in the FIR, and 
found that 972 of them contained identifiable cores, as de- 
fined by peaks in their FIR emission at 250 /im. We classified 
the cores depending on the presence or absence of an embed- 



ded MIR source at 8- or 24-/im. We identified the youngest, 
least-evolved cores as those which contained no embedded 
MIR point sources, and labelled them starless. We found 170 
IRDCs containing starless cores. 

Previous study has shown that high mass star forma- 
tion is unlikely to h ave already occurred in starless cores 
l| Wilcock et al.|[201ll ). However, it should be noted that, if 
the column density were high enough, it is possible for an 
IRDC to absorb the MIR emission from more evolved cores 
and ap pear starless despite the pr esence of embedded pro- 
tostars l|Pavlvuchenkov et alJlioill ). 

In this paper we identify the most isolated of these 
170 starless cores, based on their appearance at FIR wave- 
lengths. We find 20 starless cores that satisfy our conditions 
for being the most isolated. We then model these 20 cores 
to determine their physical parameters. 



2 DATA 

The Hersche^ Space Observatory was launched in May 
2009. Herschel carries three science instruments. These 
are: the Spectral and Photometric Imaging Receiver 



1 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with participation from NASA. 
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Figure 1. Example SEDs for cores 307.495+0.660, 318.573+0.642 and 322.334+0.561. The points show the observed fluxes at each of 
our wavelengths and the model SED is shown as a dashed line. The 70-/^m data provide upper-limits only. 



(SPIRE. [Griffin et alj|201oh : the Photodetector Array Cam- 
era and Spectrometer fPACS. IPoglitsch et al.l [2010): and 
the Heterodyne Instr ument for the Far Infrared (HIFI, 
Ide Graauw et alj|2010h . Herschel is capable of observing in 
the FIR between 55 and 671 /im. We only use data here 
from SPIRE and PACS. 

The data used in this paper were taken as part of 
the Herschel Infrared Galactic Plane Survey (Hi-GAL), an 
Open Time Key Proje ct of the Herschel Space Observatory 
ijMolinari et alj|2010al lbh. Hi-GAL aims to perform a survey 
of the Galactic Plane using the PACS and SPIRE instru- 
ments. The two are used in parallel mode to map the Milky 
Way Galaxy simultaneously at five wavelengths (70, 160, 
250, 350 and 500 ftm), with resolution up to 5" at 70 pm. 

PACS data reduction was performed usin g the Her - 
schel Interactive Pipeline Environ ment (HIPE; lot j 2010), 



with some additions described by IPoglitsch et al 



2010) 



The standard deglitching and cross-talk correction were not 
used and custom proce dures were written for drift removal 
ijTraficante et alj|201ll ). 

SPIR E data processing used the standard processing 
methods ([Griffin et al.l 120101) . with both standard deglitch- 
ing and drift removal. In both ca ses, the ROMAGAL G ener- 
alised Least Squares algorithm l|Traficante et all 1201 ll ) was 
used to produce the final maps. A more d etailed discussion 
of the data reduction process is given by iTraficante et al.l 
<|201ll ). 

We also make use of 8-/xm and 24- /im data from Spitzer^ 
that were taken as part of the Galactic Lega cy Infrared Mid- 
Plane Survey Extraor dinaire (GLIMPSE ; iBeniamin et al.l 
120031 ) and MIPSGAL (|Carev et al.ll2009l ). respectively. We 
used the mosaics available from the Spitzer Science Centre 
to create 8- and 24-/nm maps of the region 300°^Z^330° 
and |b|«Sl°. 



2 The Spitzer Space Telescope was operated by the Jet Propul- 
sion Laboratory at the California Institute of Technology under 
a contract with NASA. 



3 SOURCE SELECTION AND MODELLING 
3.1 Source sample 

We took as our in itial sample the 170 IRDCs with starless 
cores identified bv lWilcock et all ([2012) . A Gaussian profile 
was fitted to the 250-/im intensity of each core. This was 
used to determine the full-width half-maximum (FWHM) 
of the core. Each starless core was also studied at 160 and 
500 to determine if it was isolated. Our model is best 
suited to isolated cores because isolated cores are more likely 
to have the simple internal structures. 

For a core to be deemed isolated, it has to be unconfused 
by a neighbouring core within the IRDC. We defined this by 
saying that there must be a minimum depth of magnitude 
greater than 3 a between the peak of one core and that of any 
neighbouring core. This must be true at all FIR wavelengths 
(160, 250, 350 and 500 /im). Of the 170 IRDCs with starless 
cores, only 20 cores were found to be truly isolated according 
to our rigorous definition. These cores are listed in Table [1] 



3.2 Distances 

The distances to our infrared dark cores had not been pre- 
viously calculated. The galacto-ce ntric distance to each core 
was therefore estimated using the lBrand fc Blitz! l|l993h ve- 
locity curve: 



Q 1 



dn 



g,0 



+ a 3 



(1) 



where ai = 1.0069, a 2 = 0.0363 and a 3 = 0.0065 in the 
fourth quadrant of the Galactic Plane. 6 is the core's or- 
bital velocity and d g is its galacto-centric distance. #o and 
d g fi are the solar equivalents. Canonical values of 8.5 kpc 
and 220 kms -1 were used for ri Q ,p and dp re spectively. CO 
velocities were taken from iDame et all ll 19871). This method 
is only accurate up to 8 kpc ( Brand fc Blit jri993l ). 

From this, the heliocentric distance was calculated using 

d g — (d\ cos 2 b + dg — 2 d g ,o dh cos(b) cos(Z)) 3 (2) 

where dh is the core's heliocentric distance and b and I are 
the Galactic latitude and longitude of the core. This is solved 



© 0000 RAS, MNRAS 000, 000-000 



4 Wilcock, L. A. et al. 





1 





i 


i 




1.1 
1 




1 1 

+ 


i i 

+ 


i 

+- 


- 


+ 




t 

_1_ 


++ 




0.9 
0.8- 




+ 

+ 


++ 




+ 


+ 


+ 

++ 

+ 


+ 

+ 




Radius (pc) 


0.7 
0.6 
0.5 




+ 

+ 


y 




+ 


+ 










0.4 


■ + 












+ 








0.3 




+ y" 


















0.2 











3000 3500 

Distance (pc) 



3000 3500 

Distance (pc) 



Figure 2. Mass (left) and radius (right) versus core helio-centric distance. The dashed line shows the cutoff point for a core's radius 
a core with a radius below this line will not have been selected for modelling due to our 'isolated' criterion (see text). 



as a quadratic equation and results in two possible distances 
for each core — one on the near side of the Galaxy and one 
on the far side. 

The near/far distance ambiguity exists in all cases. 
However, each of our cores had to show up in absorption 
at 8 )im to be considered a candidate IRDC in the first 
place (PF09). Therefore, most of the MIR-emitting mate- 
rial must be behind our cores, rather than in front of them, 
or absorption would not be evident. For this reason, we as- 
sume that the near distance is correct for all of the cores. 
The calculated distance to each core is listed in Table [1] 
Fourth quadr ant IRDCs have a t ypical galacto-centric 



distance of 6kpc (I Jackson et al.ll2008|) and a ty pical helio- 
centric distance of 4kpc ( Peretto fc Fullei1l2010t ). The mean 



heliocentric distance to our isolated cores is 3.1 kpc, closer 
than 4 kpc. This is due to a selection bias in our isolation 
criterion. It can be seen in Figure [2] where we compare the 
masses (see discussion below) and radii of the cores to their 
distance. Only cores with a high mass or large radius were 
seen as isolated at greater distances. This is due to the an- 
nular resolution of Herschel. Using our lowest data-points, 
a cutoff is extrapolated This is shown as a dashed line in 
Figure [2] Below this cutoff a core is unlikely to be selected 
for modelling due to our isolation criterion. 



3.3 Radiative Transfer Modelling 

Our 20 isolated , starless cores were modell e d us- 
ing Phaethon (IStamatellos fc Whitworthl I20031 . 120051 ; 
IStamatellos et al.ll201fj| ). Phaethon is a 3D Monte Carlo 
radiative transfer code. The code uses luminosity packets to 
represent the ambient radiation field in the system. These 
packets are injected into the system where they interact (are 
absorbed, re-emitted or scattered) stochastically. 

The input variables of the code are the density profile, 
the strength of the interstellar radiation field (ISRF), the 
dust properties of the system, the size of the core and its 
geometry (i.e. spherical, flattened or cometary). The code 
calculates the temperature profile of the system as well as 
spectral energy distributions (SEDs - see Figure [TJ and in- 
tensity maps, at different wavelengths and viewing angles — 
see Appendix [XI 



Phaethon use s lOssenkopf fc Henningl (Il994l ) opaci- 
ties for a standard iMathis et alj ( 19771 ) grain mixture of 
52 per cent silicate and 47 per cent graphite with grains that 
have coagulated and accreted thin ice mantles over 10 years 
at densities of 10 6 c m" 3 . Further details of the dust mode l 
use d can be found inlstamatellos fc Whitworthl (|2003l. |2005[ ) 
and IStamatellos et al.l (|2010f ) . 

All the cores showed some measure of eccentricity in the 
obser vations and so were mode l led w ith a flattened geometry 
— seelstamatcll os et al. I l|2004l2010l ) for details. In this case 
the density profile is given by: 



n(r,9) = n (H 2 ) ■ 



1 + A(^) 2 [sin(0)r 



(3) 



where r is the radial distance, 9 is the polar angle and 7?o 
is the flattening radius (i.e. the radial distance for which 
the central density is approximately constant) and was set 
to one tenth of the core's semi-major axis, Rmax- ^(Hb) is 
the central density, which is controlled as an input variable. 
A is the asymmetry factor and controls the equatorial to 
polar optical depth ratio and determines how flattened the 
core is (i.e. its eccentricity), p determines how quickly the 
optical depth changes from equator to pole, and was set to 
2. These values for Ro and p were chosen as they have been 
shown previou sly to give sensible results when modellin g 
cold cores (e.g. IStamatellos et al.ll2010l ; IWilcock et al.ll201ll ). 
By restricting these parameters to sensible estimates a single 
model may be fitted to each core. 

The size and shape of each modelled core was set to 
match the observed core at 250 jj,m. Rmax was set equal to 
the FWHM of the major axis of each observed core. The 
Rmax of each core is listed in Table Q] The FWHM of the 
cores have an error of ±15 per cent, due primarily to uncer- 
tainties in the background levels. Contours from the 250 ^tm 
observations were overlaid onto images of the modelled core. 
The asymmetry factor of each model was manually varied 
until the modelled shape visually matched that of the ob- 
served core. 

To test the goodness of fit of the modelled morphology 
to the observed morphology we compared the eccentricity 
and the FWHM. For the eccentricity, a percentage differ- 
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Figure 3. Histograms of (top row, left to right) core radius, core mass, (middle row) the ISRF surrounding the core and the core's central 
density and (lower row) the surface and central temperatures found in the cores. Note that the x-axes for mass and central density use 
a logarithmic scale. 



ence of less than 20 per cent was considered a good fit. The 
eccentricities of the observed and modelled cores are given 
in Table The FWHM of the modelled core is expected 
to agree with the FWHM of the observed core within the 
prescribed errors of ~15 per cent. The FWHM of the ob- 
served and modelled cores at 250 fj,m are given in Table [3] 
Differences between the model and observed morphology are 
mainly due to the assumption that the observed cores were 
perfect ellipses and the difficulties of fitting such a model 
by-eye. 

The flux density, integrated over twice the FWHM of 
each core using an elliptical aperture, was measured in all 
five FIR maps and an SED was plotted. These flux densities 



have been background subtracted, where the background 
was defined using an off-cloud, elliptical aperture. 

The central density, no(H2), and the ISRF incident on 
the cores were manually varied until the output model's SED 
matched the observed data. Th e ISRF is tak en to be a multi- 
ple of a modified version of the lBlackl £l994) radiation field, 
which gives a good approximation to the radiation field in 
the solar neighbourhood. The final values for no(H2) and 
the incident ISRF have an uncertainty of ±15 per cent. This 
is based on a by-eye estimation of how much they can be 
varied until the model SED no longer fits the observed flux 
densities. Examples of the model SEDs are shown in Fig- 
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ure[T]as a dotted line. The final values for no(H2) and the 
incident ISRF are listed in Table Q] 

To test the goodness of fit of the modelled SED to the 
observed SED the \ 2 value was calculated for each core. 
Only the flux densities from 160-500 fim were taken into ac- 
count when calculating \ 2 ■ The 70 fim point was only used 
as an upper limit and was ignored in this calculation. Rmax 
and the asymmetry factor were set through visual obser- 
vations. As such, only the ISRF and central density were 
varied when matching the model SED to the observed SED. 
There are four data points in each SED and two free param- 
eters, meaning there are two degrees of freedom in our SED 
fit. We aim for a 5 per cent significance level and perform a 
two-sided test. For a 5 per cent significance level, with two 
degrees of freedom, \ must be between 0.051 and 7.378 (e.g. 
iRouncefield fc Holmes! ll993T ) for the model to be deemed a 
'good' fit to the data. The final \ 2 values for each model are 
listed in Table □ 

For example, we can look at the core in 307.495+0.660, 
the SED of which is shown in Figure [T] and images of the 
observed and modelled core are shown in Appendix [A] The 
model created for the core in 307.495+0.660 shows absorp- 
tion at 8 and 70 rim, comparable to the absorption seen in 
the observations. The size and shape of the model fits with 
the contours from the observed data at all wavelengths (see 
Tables [2] and [3} . The model does, however, fit slightly bet- 
ter at longer wavelengths where the observed core appears 
as a more typical ellipse. The observed FWHM of the core 
is 55 "at 160 /xm and 67 "at 500 /jm. The model core has a 
FWHM of 65 "at 160 rim and 62 "at 500 /im, in agreement 
with the observed values. The eccentricity of the modelled 
cores shows only a 13 % difference from the eccentricity of 
the observed core at 250 /im. The modelled SED fits well 
with the observed SED at all wavelengths and has a \ 2 °f 
1.41, well within our boundaries for a 5% significance level. 
The model is in agreement with the observed flux densities 
at 160, 250, 350 and 500 fim and below the upper limit set 
at 70 /im. 

Nineteen of the twenty modelled cores met the criteria 
for a good fit to the morphology and SED of the observed 
cores. The core in 326.632+0.951, however, did not meet 
our criteria for a good fit. This core has a \ 2 °f 10.21 and 
the morphology of the model showed a 50 per cent difference 
from the observed core at 250 /im. It was, nonetheless, the 
best fit possible with Phaethon. This may be indicating 
that this core has a more complex morphology than the 
model allows. 

It should be noted that the internal structure of cores 
within infrared dark clouds have not yet been observed in 
detail. Hence, our assumption of an elliptical geometry may 
be an over-simplification. If there is structure on smaller 
scales than we can resolve, this might affect our results. For 
example, small scale fragmentation might allow the ISRF to 
penetrate further into the cores, meaning that our calculated 
ISRF values may be upper limits. 

The figures in Appendix[X]show the output of the model 
at wavelengths corresponding to those of the observed data. 
The modelled images have pixels of 0.02x0.02 pc in size and 
show an area 2.56x2.56 pc in total. All of the model im- 
ages have been convolved with the appropriate telescope 
beam size : 2 "at 8 /im; 5.2 "at 70 (im; 11.5 "at 160 /mi; 18 "at 
250 /mi; 25 "at 350 /im; and 36 "at 500 /im jWerner et all 



Table 2. The eccentricities of the cores (to 1 decimal place). Col- 
umn 2 gives the eccentricity when measured from the data and 
column 3 shows the eccentricity as measured from the model. 
Column 4 is the difference between the data and model eccentric- 
ities as a percentage of the observed eccentricity to 2 significant 
figures. 



IRDC Name 


Eccentricity 


Percentage 




Observed 


Model 


Difference (%) 


305.798—0.097 


0.9 


0.9 


0.0 


307.495+0.660 


0.8 


0.9 


13 


309.079—0.208 


U.o 


u.y 


13 


309.111-0.298 


0.9 


0.9 


0.0 


310.297+0.705 


0.5 


0.6 


20 


314.701+0.183 


0.8 


0.9 


13 


318.573+0.642 


0.8 


0.9 


13 


318.802+0.416 


0.8 


0.9 


13 


318.916-0.284 


0.8 


0.9 


13 


321.678+0.965 


0.9 


0.9 


0.0 


321.753+0.669 


0.9 


0.9 


0.0 


322.334+0.561 


0.8 


0.9 


13 


322.666-0.588 


0.1 


0.1 


0.0 


322.914+0.321 


0.6 


0.5 


17 


326.495+0.581 


0.9 


0.9 


0.0 


326.620-0.143 


0.9 


0.9 


0.0 


326.632+0.951 


0.6 


0.9 


50 


326.811+0.656 


0.7 


0.7 


0.0 


328.432-0.522 


0.8 


0.9 


13 


329.403-0.736 


0.6 


0.6 


0.0 



Table 3. The FWHM of the cores at 250 /mi. Column 2 gives 
the FWHM when measured from the data and column 3 shows 
the FWHM as measured from the model. Both have errors of 
±15 per cent. See text for details. 



IRDC Name 


FWHM at 250 /im (") 




Observed 


Model 


305.798-0.097 


62 


82 


307.495+0.660 


55 


59 


309.079-0.208 


57 


65 


309.111-0.298 


38 


43 


310.297+0.705 


45 


38 


314.701+0.183 


40 


45 


318.573+0.642 


40 


44 


318.802+0.416 


45 


45 


318.916-0.284 


38 


42 


321.678+0.965 


65 


72 


321.753+0.669 


85 


79 


322.334+0.561 


55 


63 


322.666-0.588 


42 


40 


322.914+0.321 


38 


42 


326.495+0.581 


33 


39 


326.620-0.143 


73 


73 


326.632+0.951 


24 


31 


326.811+0.656 


43 


43 


328.432-0.522 


58 


57 


329.403-0.736 


47 


51 



© 0000 RAS, MNRAS 000, 000-000 



Starless cores in IRDCs 7 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 10 3 10 4 

Radius (pc) Central Density (cm" 3 ) 

Figure 4. Core mass versus radius (left) and central density (right). No correlation is seen in either case. 




Figure 5. Temperature profiles of two typical cores. The solid 
line shows core 310.297+0.705 and the dashed line shows core 
321.753+0.669. The top line of each pair shows a direction of 
8=0°, perpendicular to the midplane of the core, and the lower 
line shows 6=90 °, along the midplane of the core. 



120041 ; iGriffin et al]|201oh . For wavelengths where no emis- 
sion is visible in the model, an image showing background 
radiation is shown. Note that no attempt was made to cor- 
rectly model the surrounding area, so these images are not 
an accurate representation of the parent cloud but do show 
the cores in absorption a s see n in the observations — see also 
IStamatellos et al l (|201Ch and lWilcock et all (|201ll ). 



4 RESULTS 
4.1 Masses 

The core masses were calculated within the Phaethon 
code using an opacit y of ksooa™— 0.03cm 2 g _1 
jOssenkopf fc Hennind[l994l ). Our 20 cores have masses in 
the range ~ 7-200 Mq , with a mean of 58 Mq and a median 
of 5 2 Mq. Their distr i bution is shown in Figure [3] 

iRathborne et al.l (|2006l ) calculated the masses of 190 
infrared dark cores and found a range between 10 and 
2100 Mq, with a median of 120 Mq. However, they also 
find that 67 per cent of their cores have a mass between 30 




Surface Temperature (K) 




Central Temperature (K) 



Figure 6. Upper: The ISRF in multiples of the lBlacld l|l994l l ra- 
diation field plotted against the surface temperature. The dashed 
line shows the relation between the surface temperature and ISRF 
as described in Equation [4] Lower: The ISRF in multiples of the 
iBlaclj dl994l) radiation field plotted against the central tempera- 
ture. 



and 300 Mq. By choosing only a few of the youngest, most 
nearby cores, we are more likely to only pick those with lower 
masses, so this discrepancy is most likely a selection effect. 
Figure [4] shows the masses of our cores plotted against their 
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Mass (M Q ) 




Central Density (cm ) 



Figure 7. Upper: The mass of the cores plotted against their cen- 
tral temperatures. Lower: The central density of the cores plotted 
against their central temperatures. 



radii and central densities. The masses of the cores do not 
appear to be correlated with either parameter. 



4.2 Temperatures 

When selecting our sources the aim was to pick the youngest 
subset of infrared dark cores: those which show no inter- 
nal MIR source and which were thus deemed to be starless. 
The lack of a MIR source makes it unlikely that these cores 
have any significant heating from within. The cores are, in- 
stead, heated by external radiation impinging onto their sur- 
faces. This external heating results in a temperature gradi- 
ent throughout each core, decreasing from surface to centre. 
Temperatures i n individual cores have been observed to v ary 
by up to 15 K (|Peretto et al.ll2010l ; Iwilcock et al.ll201ll ). It 
should be noted that the temperature profiles of the cores 
are not linear but rather show an exponential relationship 
with distance from the centre of the core. Two example tem- 
perature profiles are shown in Figure 

In Table[T] the temperature at the centre of the core and 
at the surface are given. These have averages of 11 and 21 K 
respectively. The dust temperatures in flattened cores show 
a variation at different directions within the cores. The cores 
are colder in their midplanes (8 — 90°) than in the direction 
perpendicular to their midplanes (6 = 0°). The outer layers 




20 22 24 

Surface Temperature (K) 

Figure 8. The temperature at the centre of the core versus the 
temperature at the surface of the core. 



of the cores show the biggest variation with direction. The 
surface temperatures show a greater range (varying by over 
10 K) than the central temperatures (which have a range of 
less than 7K). Figure [3] shows the distribution of both sets 
of temperatures. 



4-2.1 Surface Temperatures 

Figure [6] shows the ISRF surrounding each core plotted 
against both the surface (upper panel; on a log scale) and 
central (lower panel) temperatures. The surface temperature 
is seen to be strongly dependent on the amount of radiation 
falling on the core. As a core's major source of heating is 
from external radiation, this is as expected. The higher the 
local ISRF, the more radiation will fall on the core and the 
higher its surface temperature will be. The two show an ex- 
ponential relationship, with increasing amounts of radiation 
needed to heat the core to greater temperatures. 

The correlation coefficient between the ISRF and T s is 
r=0.99, p<0.05. r is the Pearson-product moment correla- 
tion coefficient and p is the level of significance (i.e. p<0.05 
implies a 95 per cent confidence that a relationship exists) . 

The dotted line in Figure [6] shows the best-fit canonical 
relationship: 



In (ISRF) = 5.8 ± 0.2 In (T s ) - 16.8 ± O.i 



(4) 



where ISRF is the radia tion field fo und in the vicinity of the 
core in multiples of the iBlackl l| 19941 ) radiation field and T s 
is the core's surface temperature. 

From Equation [4] we extrapolate that for a core to have 
a surface temperature over 30 K the ISRF would need to be 
over 20 times that of the solar neighbourhood and for a 
surface temperature over 40 K the ISRF would have to be 
100 times greater. Most infrared dark cores will therefore 
have a surface temperature lower than 30 K. The majority 
of our cores have surface temperatures less than 25 K. 

Similarly, a core with a surface temperature less than 
10 K would need to be situated in a radiation field of only 
one-hundredth of the solar value. This makes lower tempera- 
tures very unlikely, although not impossible. All of our cores 
had surface temperatures greater than 16 K. 
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Figure 9. Upper: The central density of the cores plotted against 
the difference between the surface and cen tral te mperatures (AT). 
Middle: The ISRF in multiples of the iBlacld lll994D radiation 
field plotted against AT. Lower: The masses of the cores plot- 
ted against AT. 



4-2.2 Central Temperatures 

The temperatures at the centres of the cores vary according 
to different criteria to those at the surface. The central tem- 
peratures of the cores do not show any dependency on the 
ISRF - see the lower panel of Figure[6] Instead, the temper- 
ature at the centre of the core appears to be determined by 
its mass and central density, as shown in Figure [7] 

The central temperature shows an inverse linear rela- 



tionship with central density with a correlation coefficient 
of r=-0.46, p<0.05. As radiation falls onto the core surface 
it will heat the outer layers first, giving the temperature 
gradient visible in all IRDC cores. In denser cores the radia- 
tion will have a greater probability of being absorbed in the 
outer layers and less will travel through the core to heat the 
centre. Thus, denser cores have a lower central temperature 
than cores with a lower density. 

The masses of the cores show no relation to the tem- 
perature at the surface of the core but do show a weak anti- 
correlation with the central temperatures, r=-0.35, p<0.05. 
Cores with high masses have a lower central temperature. 
This is due to the same effect that causes the central temper- 
ature to vary with density. As mentioned in Section fOI mass 
and central density show no correlation with one another. 
This implies that their relationships with central tempera- 
ture must be independent of each other. 

Figure [8] shows the central temperature plotted against 
the surface temperature. The two show a broad correla- 
tion (r=0.44, p<0.05), implying that cores with a similar 
central density, have a higher central temperature if there 
is a greater ISRF. However, the cores with highest sur- 
face temperatures do not necessarily have the highest cen- 
tral temperatures. For example, cores 307.495+0.660 and 
326.811+0.656 have vastly different surface temperatures 
but identical central temperatures. 



4-2.3 Temperature Range 

The range of temperatures in each core, AT, is found via 

AT=T S -T C , (5) 

where T c is the central temperature of the core, and T s 
is the surface temperature of the core. AT is indicative of 
a core's temperature profile. A core with a high AT will 
show a steeper profile than a core with a lower AT. Figure 
[5] shows the temperature profiles of two cores with similar 
radii. G310. 297+0. 705 has AT= 10.5 K and a far steeper 
profile than G321. 753+0. 669 which has AT= 4.8 K. AT ap- 
pears to be dependent on a core's central density (r=0.81, 
p<0.05) and the surrounding ISRF, but not with its mass 
(r=0.27, p>0.05). This is shown in Figure [§] A possible ex- 
ponential relation between AT and the ISRF can be seen. 
As this is similar to the relationship between ISRF and the 
surface temperature, it is likely a result of the AT depen- 
dence on surface temperature of the core. 

The strongest relation is that between AT and the core 
central density. The difference between the temperature at 
the surface and centre of a core clearly increases as the cen- 
tral density increases. Radiation passing through denser ma- 
terial has a greater probability of being absorbed in the sur- 
face layers of the core, and so less radiation passes through 
to (and heats) the centre of the core. Therefore, if two cores 
are situated within areas of equivalent ISRFs then, while 
the surface temperatures would be the same, the denser core 
would have a colder centre and exhibit a wider range of tem- 
perature. 

AT also shows a weak inverse correlation with the ra- 
dius of a core (r=-0.55, p<0.05) - see Figure [TO] The larger 
the size of the core, the smaller the temperature difference 
between the surface and the centre. This is the opposite of 
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Figure 10. AT versus radius. 
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Figure 11. Central density versus radius. 



what one would generally expect and, as it is a weak rela- 
tionship, is not likely to be a direct result but rather a side 
effect of the fact that our smallest cores have the highest 
densities. 



4.3 Central Densities 

The central densities of the cores vary from 1.5xl0 3 cm -3 
to 2.0xl0 5 cm~ 3 , with a mean of 2.0xl0 4 cm~ 3 . The cen- 
tral density of the core is an input to the Phaethon code 
and mostly affects the longer wavelength end of the SED. 
Figure[3]shows the distribution of central densities. The den- 
sity profile of each core is calculated via Equation [3] This 
results in the density at the surface of the core (where it is 
at its lowest) directly correlating with the maximum, central 
density. It is simply an input definition of the model. The 
surface densities range from 0.1 to 20cm -3 . 

Figure [TT1 shows the central density of each core plotted 
against its radius. There is no obvious correlation between 
the two but there is a noticeable lack of cores in the lower 
left-hand corner with both a low central density and a small 
radius. This is likely to be a selection effect. This implies 
that: either small, low density cores do not show up as a 
distinct peak against their parent cloud; or they do not block 
enough MIR background emission to show a significant dip 



at 8 /j,m; or that they are simply fainter and so are more 
likely to be confused at 500 /im. 

This latter theory was tested by creating two model 
cores that would reside in the lower left of Figure 1111 Both 
had a radius of 0.3 pc, equal to the smallest core we model. 
One had a central density of 2 x 10 3 cm -3 and the other had 
2 x 10 4 cm -3 . We placed them both at a di s tance of 3.1 kpc, 
gave them an ISRF of 3.2 times the IBlackl (|l994l ) radiation 
field and an asymmetry factor of 2.5 (these are the average 
parameters of our modelled cores). The resulting objects 
have masses of 0.4 Mq and 3.9 Mq. The cores were added to 
the Hi-GAL data in several regions and studied as with the 
original PF09 cores. Even in very low background regions, 
neither core met our criteria for isolation at 500 /im. They 
would not have been selected by our observations. 

4.4 The Interstellar Radiation Field 

The ISRF fo r each core is quote d in terms of m ultiples of the 
IBlackl (| 19941 ) radiation field. The lBlackl (|l994T l radiation field 
models the ISRF in the solar neighbourhood and measures a 
value of 1.3 x 10" 16 erg s" 1 cm -2 Hz -1 sr" 1 at 500 (im. Along 
with the central density, the ISRF is one of the inputs of the 
Phaethon code (see Section I3.3|l , and mostly affects the 
SED at shorter wavelengths. T he ISRF ar ound our cores 
ranges from 0.6 to 10.4 times the lBlackl (|l994T ) field, with an 
average multiple of 3.2. Figure [3] shows the distribution of 
the cores' ISRF. 

The amount of radiation incident on our cores is typ- 
ically greater than that which we observe in the solar 
neighbourhood. This is consistent with the inner Galactic 
Plane generally having a stronger radiation field than that 
of the solar neighbourhood. However, two of the cores do 
show a radiation field lower than the solar value. These are 
307.495+0.660 and 321.678+0.965 (see Table There are 
two possible explanations for this: either they simply exist 
in an area which has a particularly low ISRF; or the ISRF 
surrounding the parent IRDC is consistent with other cores 
but the parent cloud is absorbing most of the radiation. The 
latter could be caused by the IRDC being especially dense, 
or by the cores being more deeply embedded within the par- 
ent cloud. Either case would make it difficult for radiation 
to penetrate far enough to affect the cores. 

Figure [121 shows a plot of core radius versus ISRF. We 
detect cores of all sizes with ISRFs at the lower end of our 
scale but our modelled sample doesn't include any cores with 
a high ISRF and a large radius. If radiation pressure were 
forcing the cores to condense then we would expect the cen- 
tral density to increase as the radius decreased. While our 
smaller cores do tend to have densities at the higher end of 
our range, the larger cores have both low and high densities. 
This is due to the selection effect discussed in Section 14.31 
We therefore dismiss radiation pressure as the reason why 
we detect no large cores in areas of high ISRF. Instead, we 
consider that in areas with a high radiation field the outer 
layers of large cores may 'boil away' due to the large amounts 
of energy being absorbed. This would explain why we only 
see small cores in areas of high ISRF and why there is no 
significant increase in the density of these cores. 

The radius and surface temperature appear to show a 
similar, albeit inverse, relationship to that of ISRF and ra- 
dius. We speculate that the radius of a core is actually af- 
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Figur e 12. Core radius versus ISRF - in multiples of the iBIackl 
( 1994) radiation field. 



Table 4. The mean values of the physical properties of the 20 
modelled cores. 



Physical Property 




Mean 


Radius (pc) 




0.7 


Central Density (cm" 




5.6xl0 4 


Mass (Mq) 




58.3 


ISRF 




3.2 


Central Temperature 


(K) 


11.5 


Surface Temperature 


(K) 


21.0 


AT (K) 




9.5 


Asymmetry Factor 




2.5 



fected by the ISRF and that any correlation between tem- 
perature and radius exists as a result of the dependence of 
the surface temperature on the ISRF. 



5 SUMMARY 

Using data from Herschel and Spitzer, 20 isolated, starless 
cores were modelled with the Phaethon radiative transfer 
code. A flattened density profile, as described by Equation^ 
was assumed for all. The FWHM of each core was measured 
at 250 fim and used as the core's radius. The radial distance 
for which the central density is constant was set to one tenth 
of this value. The density and surrounding ISRF were varied 
until the model SED matched the observed SED at 160, 250, 
350 and 500 /im. 

Output parameters were measured from the model. The 
IRDC cores were found to have masses ranging from around 
ten to around a few hundred solar masses. The masses of the 
cores were found not to correlate with their size or central 
density. The temperature at the surface of a core was seen to 
depend almost entirely on the level of the ISRF surrounding 
the core. No correlation was found between the tempera- 
ture at the centre of a core and its local ISRF. The central 
temperature was seen to depend, instead, on the density and 
mass of the core. A core with a high density or mass is likely 
to have a lower central temperature. 

The range of temperatures from centre to edge in any 
one core was found to correlate strongly with the central 
density of that core. Low density cores exhibit a much 
smaller range of temperatures and thus a shallower temper- 
ature gradient than higher density cores. The cores situated 
in a high ISRF (and thus with a high surface temperature) 
are also likely to show a greater range of temperature from 
centre to edge. 

No large cores were found within an area of compara- 
tively high ISRF. However, there appears to be no correla- 
tion between the ISRF surrounding a core and its central 
density. Hence, we rule out radiation pressure as the cause. 
Instead, we consider that the outer layers of cores in areas 
of high radiation may be evaporated due to the high levels 
of energy being absorbed, leaving only smaller cores in high- 
ISRF regions. 16 out of 20 isolated IRDC cores showed high 
levels of ellipticity. 



4.5 Asymmetry Factor 

The asymmetry factor input to the model controls the shape 
of the core. A value of 1.0 would result in a circular core. The 
ellipticity of the core is increased by raising the asymmetry 
factor up to a maximum of 3.0. Most IRDC cores are elon- 
gated and have a high eccentricity. 80 per cent of our cores 
have an asymmetry factor greater than 2 - see Table [T] It is 
also possible that some of our more circular cores could be 
elongated but viewed end-on. 



4.6 Mean Core Properties 

Table [3] shows the mean values of each of the eight physical 
properties discussed in Sections 14. 1H4. 51 The mean physical 
properties reiterate the fact that starless cores are dense, 
elongated, cold objects and, when found in the Galactic 
Plane, have a local ISRF greater than that which is found in 
the solar neighbourhood. Within the cores, the temperature 
decreases from surface to centre, typically by ~10K. The 
mean asymmetry factor is 2.5, corresponding to an aspect 
ratio of ~1:7. 
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Starless cores in IRDCs 13 




Figure Al. 305.798-0.097 - Upper row: Observational images of the infrared dark core at (left-right) Spitzer 8 /im, PACS 70 fim and 
PACS 160 fim. Second row: Modelled images at 8 /im, 70 fim and 160 (im. Third row: Observational at SPIRE 250, 350 and 500 fim. 
Lower row: Modelled images at 250, 350 and 500 fim. All images are ~0. 05x0. 05° in size, equivalent to 2.56x2.56 pc. 
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